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ABSTRACT 



X 



Of the many ways of detecting high redshift galaxies, the selection of objects due to their redshifted Lya emission has become one 
of the most successful. But what types of galaxies are selected in this way? Until recently, Lya emitters were understood to be small 
star-forming galaxies, possible building-blocks of larger galaxies. But with increased number of observations of Lya emitters at 
lower redshifts, a new picture emerges. Lya emitters display strong evolution in their properties from higher to lower redshift. It has 
previously been shown that the fraction of ultra-luminous infrared galaxies (ULIRGs) among the Lya emitters increases dramatically 
between redshift three and two. Here, the fraction of AGN among the LAEs is shown to follow a similar evolutionary path. We argue 
that Lya emitters are not a homogeneous class of objects, and that the objects selected with this method reflect the general star forming 
and active galaxy populations at that redshift. Lya emitters should hence be excellent tracers of galaxy evolution in future simulations 
and modeling. 
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1. Introduction 

Several decades of studies of high redshift galaxies have shown 
that the star formation rate density of the Universe peaked 
around redshift z ~ 2 (e.g. Hogg et al. 1998, Hopkins 2004, 
Hopkins & Beacom 2006). At the peak of the star formation his- 
tory, nearly ten times more stars were formed than in our local 
Universe, whereas at higher redshifts, the star formation density 
was equally low as it is now. Similarly, a trend in the volume 
density of AGN has been found with a peak at z = 1.5 — 2 (e.g. 
Miyaji et al. 2000, Wolf et al. 2003, Bongiorno et al. 2007). The 
coincidence that the two density functions, for star formation and 
numbers of AGN, peak at similar redshifts has been proposed as 
being due to both of these properties being linked to the hier- 
archical build-up of galaxies and mergers of dark matter haloes 
(e.g. Kauffmann & Haehnelt 2000, Bower et al. 2006). 

As for the high redshift Universe, one of the strongest emis- 
sion lines observable is the Lyman-a (Lya) line. By now, several 
hundreds of Lya emitters (LAEs) have been detected through 
narrow-band imaging at z = 0.3 — 7.7 (e.g. M0ller & Warren 
1993, Fynbo et al. 2003, Gronwall et al. 2007, Venemans et 
al. 2007, Nilsson et al. 2007, Finkelstein et al. 2007, Ouchi et 
al. 2008, Grove et al. 2009, Hibon et al. 2010). Lya emission 
may be generated by three main mechanisms; the ionising flux 
of O and B stars, indicative of star formation, the ionising flux of 
an energetic UV source, e.g. an active galactic nucleus (AGN), or 
due to infall of gas on a massive dark matter halo (c.f. Dijkstra et 
al. 2006a,b, Nilsson et al. 2006). The volume density of sources 
where the Lya emission is dominated by the latter is expected to 
be very low compared to those where the Lya emission comes 
from star formation or AGN sources, hence, the volume density 
of Lya emitting objects found in the Universe is expected to fol- 
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low the general evolutionary occurrences of the star formation 
history, and the AGN history, with redshift. 

In this Letter we discuss the fractions of ULIRGs and AGN 
among Lya emitters at different redshifts. In a previous publica- 
tion, the ULIRG fraction among the LAEs was shown to exhibit 
a sharp transition from very few to a larger sub-sample at a red- 
shift around 2.5 (Nilsson & M0ller 2009). Here, the fraction of 
AGN among LAEs is shown to follow a very similar relation, 
indicating that the underlying galaxy population is transition- 
ing rapidly from z > 3 to z ~ 2, (see also a similar result in 
Bongiovanni et al. 2010). We here ask the question how these 
results relate to the general galaxy evolution in the Universe. 

Throughout this paper, we assume a cosmology with Hq = 
72 km s- 1 Mpc"\ Q m = 0.3 and n A = 0.7. 



2. The AGN fraction of LAEs 

To determine AGN fractions among LAEs at different redshifts, 
we started with the data-set of Nilsson et al. (2009, 201 1). This 
sample of LAEs was found by narrow-band imaging with the 
ES02.2m/MPG telescope on La Silla, Chile, using the Wide- 
Field Imager (WFI). A central section of the COSMOS field was 
searched for Lya emitters with z — 2.206 — 2.312. Details of the 
data reduction and candidate selection can be found in Nilsson 
et al. (2009), where in total 187 LAE candidates were found. 
Follow-up spectroscopy of 152 candidates was performed with 
the VLT/VIMOS instrument early 2010. The results of this cam- 
paign will be published in a forthcoming publication, but we will 
here exclude those candidates which were not confirmed in the 
spectroscopy. This brings the total sample of LAEs to 171. AGN 
were selected among the LAEs by use of X-ray data from XMM 
(Hasinger et al. 2007), Chandra (Elvis et al. 2009) and the VLA 
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Fig. 1. Cumulative fraction of non-AGN LAEs as a function of 
Lya luminosity. Black points are from the Nilsson et al. (2009) 
sample. Coloured points are from publications at redshifts z = 
0.3 (Atek et al. 2009, Cowie et al. 2010, Finkelstein et al. 2009c, 
Scarlata et al. 2009), z = 2.1 (Guaita et al. 2010), z = 3 
(Gronwall et al. 2007, Ouchi et al. 2008) and at z ~ 4 - 5 (Wang 
et al. 2004, Ouchi et al. 2008). Lines are best fit models (see text 
for details). 



Fig. 2. Evolution in AGN fraction above a given Lya luminosity 
as a function of age of the Universe. Colour scheme is identical 
to Fig. Q] and refers to redshift. Different symbols mark the AGN 
fractions above three representative Lya luminosities. The black 
lines are fits to the data, see text for explanation. A transition red- 
shift at z ~ 2.5 (age ^2.5 Gyrs) are seen in all fits independent 
of luminosity limit. The purple line shows the redshift transi- 
tion of the fraction of ULIRGs in LAE samples as presented in 
Nilsson & M0ller (2009). 



(Schinnerer et al. 2007). For details, see Nilsson et al. (2009, 
201 1). In total, 24 LAEs are selected as AGN. 

It was shown already in Nilsson et al. (2009) that the AGN 
fraction among these lower redshift galaxies is larger than in 
samples of LAEs at higher redshifts. With the further Chandra 
detections presented in Nilsson et al. (201 1), this surplus of AGN 
is even larger. However, since comparing AGN fractions from 
different surveys has a bias due to the Lya luminosity reached, 
a more careful investigation requires an analysis independent of 
this bias. In Fig.[T]the AGN fraction is thus presented as a func- 
tion of Lya luminosity. The black points are from the Nilsson 
et al. (2009) sample. It is clear that above some flux limit, all 
Lya emitting objects are quasars. Similarly, the function tends 
towards one at fainter fluxes, where the non-AGN LAEs will 
outnumber the AGN LAEs. 

In Fig. Q] the results of several other narrow-band surveys 
for LAEs are also shown. As we have no access to the exact 
Lya fluxes of the AGN in most samples, they are shown as sin- 
gle, cumulative points at the Lya flux limit of each survey. The 
results include AGN fractions from redshifts z = 2.1 (Guaita 
et al. 2010), z ~ 3 (Gronwall et al. 2007, Ouchi et al. 2008), 
z = 4 - 5 (Wang et al. 2004, Ouchi et al. 2008), and at low red- 
shift, z = 0.3 (Atek et al. 2009, Cowie et al. 2010, Finkelstein et 
al. 2009c, Scarlata et al. 2009). 

The data-points at z ~ 2.2 (the data from the z = 2.3 sam- 
ple, and from Guaita et al. 2010) were found to be well fitted 
with an exponential AGN function of the form: 

X AGN (log L Lya ) = exp u -— i — (1) 

\ -X-AGNfl J 

In this equation, the two free parameters are the L^^agn which 
defines the Lya luminosity where the AGN fraction goes to 
one, or the other way around that the non-AGN LAE fraction 
among the narrow-band selected sample goes to zero, and the 
normalisation constant Xagn.o that determines how quickly 



the AGN function goes to zero. A further criteria was that the 
AGN fraction is unity if log(L,L ya ) > logfX^cw), i.e. that all 
Lya emitting objects above a certain Lya luminosity are AGN- 
powered. The fits to the z ~ 2.2 sample is shown in Fig. Q] 
with best fit parameters log(L*,AGN) — 43.61*0 27 er 8 s 1 
and X A gn,q = 0.59±g;?J. For the other redshifts, the L* agn 
was kept constant and only the normalisation constant was deter- 
mined using all data-points at a particular redshift. This is based 
on the two assumptions that i) the functional form of the AGN 
fraction is similar at all redshifts, and ii) that the Lya luminosity 
above which all emitters are AGN does not change with redshift. 
These fits are also shown in the Figure. A clear evolution from 
higher to lower redshift appears to be present. 

3. Evolution in the AGN fraction with time 

To illustrate the evolution of the LAE AGN fraction with time, 
the fraction of AGN LAEs in the sample was extracted from the 
fits at three different Lya luminosities, i.e. the AGN fraction that 
a given exponential fit in Fig.[T]gives for a certain Lya luminos- 
ity, and plotted as a function of redshift in Fig. [2] Worst case 
uncertainties on these values were found by varying the expo- 
nential fits within their ler uncertainties and finding the mini- 
mum/maximum AGN fraction expected. In this Figure, a clear 
trend towards larger AGN fractions at lower redshifts is seen, 
and in particular a rather sudden transition from very small frac- 
tions at z > 3 to significant fractions at z < 2. This trend mimics 
the trend of increasing ULIRG fraction among the LAEs with 
decreasing redshift, as found in Nilsson & M0ller (2009; also 
plotted in Fig. 01. The AGN data were fitted, minimising the \ 
function, with the same transitional equation as in Nilsson & 
M0ller (2009): 

A F 

AF(z) = _£(l-tanh(0 (z - z tr ))) (2) 
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Table 1. Best fit results for Eq|2] 



log L LyQ 
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0.311°;™ 
<n+ - 09 
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i ti+0.59 
1 - :,i -0.77 
Q1+ - 67 


ULIRGs 


2.52 


0.2 


2.28 



Notes. The results at the bottom of the table refer to the best 

fit parameters for a sample of LAE ULIRGs, presented in Nilsson 
& M0ller (2009). 



In this equation, AF is the AGN fraction, 8 represents the steep- 
ness of the transition and z tr is the transition redshift. The best 
fit values for the different Lya luminosity limits can be found in 
TableQ] The most interesting parameter is the transition redshift. 
It appears that the brightest AGN transition first, while fainter 
AGN remain less numerous until later times. This is consis- 
tent with a down-sizing scenario, where the larger galaxies/AGN 
form first (Juneau et al. 2005). Comparing the results with those 
for the LAE ULIRGs shows the transition for the ULIRGs to be 
sharper than those for the AGN, with a transition redshift nearer 
to those of the brighter LAE AGN. 

4. Discussion 

4. 1. AGN and LAE volume density evolution 

Having determined the AGN and ULIRG fractions among 
LAEs, one may ask the question how these fractions relate to 
the general evolution of galaxy properties with redshift. An 
independent test of the determined AGN fraction is to cal- 
culate the LAE volume density, by dividing a known AGN 
volume density by the LAE AGN fraction. For this test, we 
use the AGN comoving volume densities as determined ob- 
servationally by Bongiorno et al. (2007) and theoretically by 
Hopkins et al. (2007), and we compare with comoving LAE 
volume densities from eighteen different publications (Steidel et 
al. 2000, Stiavelli et al. 2001, Dawson et al. 2004, Hayashino et 
al. 2004, Shimasaku et al. 2006, Nilsson et al. 2007, Gronwall et 
al. 2007, Venemans et al. 2007, Murayama et al. 2007, Prescott 
et al. 2008, Deharveng et al. 2008, Ouchi et al. 2008, Nilsson et 
al. 2009, Grove et al. 2009, Shioya et al. 2009, Guaita et al. 2010, 
Bongiovanni et al. 2010, Yuma et al. 2010). In Fig.[3]these LAE 
volume densities are shown, together with the predictions. In the 
following, all Type 1 AGN are assumed to have strong enough 
Lya emission lines to be detected as LAEs. The predictions in 
the left panel are done for AGN brighter than Mb < —22, as 
the z = 2.25 LAE AGN from Nilsson et al. (2009) are all 
brighter than this magnitude. The scatter in the data-points is 
large, probably due to different selection criteria for the differ- 
ent surveys, but the overall trends are clear; the LAE volume 
density remains relatively flat over the redshift range z = 2 — 6, 
but falls at low redshift. The volume density of LAEs predicted 
by the AGN fraction, in turn, has a steep increase towards higher 
redshift. In the range z — 2 — 3, the predictions and the mea- 
surements agree well, but there is an apparent lack/surplus of 
LAEs at high/low redshift. One explanation to this trend may be 
that the LAE AGN have different luminosities at different red- 
shifts. In the right panel of Fig. [3] the effect of varying the AGN 
luminosity is shown. It becomes clear that the lack of LAEs at 
very high redshift (z > 4) can be explained if the AGN are all 



very luminous, and similarly that the surplus of AGN at z = 0.3 
can be explained if all AGN are faint. The fact that the LAE 
volume density agrees so well with the predictions is an inde- 
pendent confirmation of the derived AGN fraction evolution, as 
well as an indication that the LAE AGN follow the general AGN 
evolutionary trends in the Universe. 

4.2. Lya emitters as tracers of galaxy evolution 

The apparent trends displayed by LAEs as a function of redshift 
are easily explained with a simple, phenomenological model 
when considering which objects in the Universe emit Lya pho- 
tons. At high redshift (z > 3), most galaxies have low stel- 
lar masses and are in their first star forming episode. AGN and 
ULIRG number densities are low. Most Lya photons at this time 
will thus trace this; young, small star-forming galaxies in their 
first starburst, with little dust content. Very few AGN or ULIRG 
LAEs will be found. A Gyr later, at z ~ 2, both the star forma- 
tion rate density and the AGN number density are peaking. Some 
galaxies are already in their second, or later, starburst. Now the 
sample of Lya selected galaxies will reflect this diversity, with 
more AGN and ULIRG LAEs detected and a larger range of 
stellar properties (cf. Pentericci et al. 2009, Nilsson et al. 201 1). 
Whereas there will be some galaxies in the same stage of evolu- 
tion as at higher redshift, there will also be some galaxies with 
evolved stellar populations that are experiencing more recent 
star formation. These galaxies can be more massive and/or more 
dusty. The study of the redshift evolution of LAEs can thus be 
seen as an evolution in Lya emission in the Universe, and could 
as such be an interesting test for galaxy evolution models includ- 
ing both star formation, dust effects and AGN number densities. 
Note that in this model, no assumption about the evolution of 
individual LAEs is made. It is unlikely that LAEs at higher red- 
shift will appear as such at lower redshift. The argument here 
is rather that LAEs are a random sub-set of all star-forming or 
active galaxies at a certain time, that happen to be in a Lya emit- 
ting phase. Also, from this data it is not possible to draw any 
conclusions on the evolution of passive (i.e. non-star-forming or 
inactive) galaxies. 



5. Conclusion 

Based on the apparent evolution of the AGN and ULIRG frac- 
tions of LAEs, fractions that are expected to follow the evolu- 
tion of star forming and active galaxy populations throughout the 
Universe, some conclusions can be drawn regarding the general 
galaxy evolution scenario. From the evolutions of the fractions 
and volume densities in Fig. [2] and [3] it is seen that at first, in the 
very young (z > 4) Universe, only stars formed, regardless of 
AGN or ULIRGs. Later, at an age of approximately 3 Gyrs or 
redshift ~ 2.5, a secondary process started, suddenly increasing 
the number of both AGN and ULIRGs in the galaxy popula- 
tion. After this abrupt increase in dusty and active galaxies, an 
equilibrium was reached that has lasted to this day. A possible 
explanation for this very stable equilibrium may be due to feed- 
back effects controlling the SFR in the galaxies. It is clear that 
this remains speculation until further data is collected, but proves 
the importance of studying objects emitting Lya emission in the 
Universe. 

Acknowledgements. The authors wish to thank Lucia Guaita for providing in- 
formation about the AGN in CDFS and Yujin Yang for useful comments on the 
draft. 
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Fig. 3. LAE volume densities as a function of redshift and limiting Lya luminosity. The data-points are collected from 18 publi- 
cations, searching for LAEs, with colour coding referring to the limiting Lya luminosity in each survey (where [blue, yellow, red] 
have log Lya [< 42.25, 42.25 - 42.75, > 42.75] erg s" 1 ). Models (see text) come from Bongiorno et al. (2007; B07) and Hopkins 
et al. (2007; H07). To the left, all measurements, and the two models (with Mb < —22.0), are shown to the different Lya luminosity 
limits. Beyond z ~ 4.5 both models and measurements become uncertain. At z — 2 — 3 the measurements agree reasonably well 
with the expectations. The volume density at z = 0.3 is significantly too large, compared to the models. To the right, only models 
and measurements to the faintest luminosity limits are shown, and instead the effect of varying the luminosity of the AGN is shown. 
For each model the gray lines show the predictions for Mb < —20.0 (above the original line) or Mb < —24.0 (below the original 
line). The volume density at z = 0.3 agrees better now, if all AGN at this redshift are very faint in the continuum. 
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